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We present a consistent way of treating a double counting problem unavoidably arising within the 
LDA+DMFT combined approach to realistic calculations of electronic structure of strongly correlated systems. 
The main obstacle here is the absence of systematic (e.g. diagrammatic) way to express LDA (local density 
approximation) contribution to exchange correlation energy appearing in the density functional theory. It 
is not clear then, which part of interaction entering DMFT (dynamical mean-field theory) is already taken 
into account through LDA calculations. Because of that, up to now there is no accepted unique expression 
for the double counting correction in LDA+DMFT. To avoid this problem we propose here the consistent 
LDA'+DMFT approach, where LDA exchange correlation contribution is explicitly excluded for correlated 
states (bands) during self-consistent band structure calculations. What is left out of Coulomb interaction for 
those strongly correlated states (bands) is its non-local part, which is not included in DMFT, and the local 
Hartree like contribution. Then the double counting correction is uniquely reduced to the local Hartree con- 
tribution. Correlations for strongly correlated states are then directly accounted for via the standard DMFT. 
We further test the consistent LDA'+DMFT scheme and compare it with conventional LDA+DMFT calculat- 
ing the electronic structure of NiO. Opposite to the conventional LDA+DMFT our consistent LDA'+DMFT 
approach unambiguously produces the insulating band structure in agreement with experiments. 
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1. INTRODUCTION 

During the last 15 years the so called LDA+DMFT 
approach (local density approximation + dynamical 
mean-field theory) became a common tool to describe 
band structure of real strongly correlated materials 
[1, 2, 3, 4, 5, 6]. In this approach the results of LDA 
band structure calculations are supplemented with local 
Coulomb (Hubbard) interaction term for those states 
which are counted as strongly correlated. Formally the 
LDA+DMFT Hamiltonian can be written as 
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Here U^ m , are the most important matrix elements of 
Coulomb matrix (Coulomb repulsion and z-componcnt 
of Hund's rule coupling) and J m m' are spin-flip terms of 
Hund's rule couplings between the strongly correlated 
electrons (assumed here to be d-states, enumerated by 
i = id and I = Id)- The prime on the sum indicates 
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that at least two of the indices of operators have to be 
different, and a =|(t) for a =t(4-). 

The LDA part of the Hamiltonian (1) is given by: 



2m e 
5p(r) ' 
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(2) 



where A is the Laplace operator, m e the electron mass, 
e the electron charge, and 



(3) 



denote the one-particle potential due to all ions i with 
charge eZi at given positions Ri, and the electron- 
electron interaction, respectively. 

The E^ A (p(r)) in (2) is a function of local charge 
density which approximates true exchange correlation 
functional £J xc [/o] of density functional theory in the 
framework of local density approximation [7] . The form 
of the function E^ A (p(r)) is usually calculated from 
perturbation theory [8] or numerical simulations [9] of 
the "jellium" model with Vi on (r) = const. Once we 
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choose some basis set of one-particle wave functions ipi 
(e.g. to do practical calculations and explicitly express 
matrix elements of the Hamiltonian (2)), we can obtain 
p as: 

N 

p(r)=5> 4 (r)| 2 . (4) 

i=i 

Finally a term H DC is subtracted in Eq. (1) to 
avoid double-counting of those contributions of the lo- 
cal Coulomb interaction already contained in i?LDA via 
Hartree term and E^ A (p(r)). Since there does not ex- 
ist a direct microscopic or diagrammatic link between 
the model (Hubbard like) Hamiltonian approach and 
LDA it is not possible to express rigorously in 

terms of U, J and p. Thus there is no unique and ac- 
cepted expression for H DC (see e.g. Ref. [10]). 

One popular expression for H DC is the Hartree like 
(fully localized limit) expression [11]: 

H DC = hjn d (n d - 1) - ij^n dff (n dCT - 1). (5) 

Here, n da = J2 m n ^dmo- = E ro (na jmff ) is total num- 
ber of electrons on interacting orbitals per spin, n d = 
J2a n dai U is Coulomb (Hubbard) repulsion and J is the 
exchange or Hund's rule coupling obtained usually from 
constrained LDA procedure [12]. The n d value can be 
obtained either from LDA calculations or can be recal- 
culated during the DMFT loop. Practically, the values 
obtained are pretty close to each other. 

Below we introduce the consistent LDA'+DMFT ap- 
proach, which allows one to avoid the double counting 
problem unambiguously. To illustrate the advantages 
of this new approach we shall apply it to calculations 
of the band structure of the well known prototype of 
charge transfer insulating system NiO. 

2. CONSISTENT LDA'+DMFT APPROACH 

One of the possible ways to solve the double count- 
ing problem is to perform Hartrcc+DMFT or Hartrec- 
Fock+DMFT calculations (see for the overview of the 
concept Ref. [13]). This approach uses the advantage 
of knowledge of diagrammatic expression for Hartree 
or Hartree-Fock terms. Thus, performing Hartree-Fock 
band structure calculations for real materials we do 
know exactly what portion of interaction is, in fact, ex- 
plicitly included. Then obviously, the double counting 
term should be chosen in the form of Eq. (5). However, 
up to now we are unaware of any Hartrcc+DMFT or 
Hartree-Fock+DMFT calculations for real materials. 

In fact, Hartree-Fock band structure calculations are 
in some sense a large step backwards from DFT/LDA 



approach, which was so successful in description of many 
real materials. Even in the case of strongly correlated 
systems DFT/LDA is recognized as a best starting point 
for further model Hamiltonian treatments, such as e.g. 
LDA+DMFT method. 

In view of this we suggest a kind of compromise be- 
tween Hartree-Fock and DFT/LDA starting points to 
be followed by DMFT calculations. As described above 
main obstacle to express double counting term exactly 
is exchange correlation E^ A (p(r)) portion of interac- 
tion within LDA. It seems somehow inconsistent to use 
it to describe correlation effects in narrow (strongly cor- 
related) bands from the very beginning, as these should 
be treated via more elaborate schemes like DMFT. To 
overcome this difficulty for these states, we propose to 
redefine charge density (4) in E^ A as follows: 

p'(r) = ^|^(r)| 2 (6) 

excluding the contribution of the density of strongly cor- 
related electrons. Then this redefined p'(r) is used to 
obtain E^ A and perform the self-consistent LDA band 
structure calculations for correlated bands at the initial 
stage of LDA+DMFT, while correlations of d-electrons 
are left to be treated via DMFT. This means that what 
is left for correlated states out of interaction on the 
LDA stage would be just the Hartree contribution of 
Eq. (2). At the same time all other states (not counted 
as strongly correlated) are to be treated with the full 
power of DFT/LDA and full p in E^ C DA . Now, the prob- 
lem of double counting correction is uniquely defined - 
it should be taken in the form of the Hartree like term, 
given by Eq. (5). 

This approach to describe realistic strongly corre- 
lated systems we shall call the consistent LDA'+DMFT. 
It is in precise correspondence with the standard defini- 
tion of correlations, as interaction corrections "above" 
Hartree-Fock. We explicitly exclude contributions to 
E^® A from (strongly) correlated bands, where correla- 
tions are treated via DMFT, while we take all electrons 
into account in LDA calculations for all other (non cor- 
related) bands. 

3. RESULTS 

Following many recent works [10, 14, 15] (and refer- 
ences therein) we choose as a testing system the proto- 
type charge transfer insulator NiO. LDA band struc- 
ture calculations for NiO were performed within the 
linearized muffin-tin orbitals (LMTO) basis set [16]. In 
the corresponding program package TB-LMTO v. 47 the 
i?x C DA was taken in von Barth-Hedin form [8]. 
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In the Fig. 1 we present LDA densities of states (left 
panel) and band dispersions (right panel) of NiO. Band 
dispersions consist of two separate sets of bands: the 
0-2p bands (from -3 to -9 eV) and Ni-3d bands, cross- 
ing the Fermi level (from 1.5 to -3 eV). Dashed lines 
in Fig. 1 show conventional LDA results. Full lines 
correspond to LDA' calculations without E^ A on Ni- 
3d states, namely, with redefined charge density (6) in 
££da Overall chan ges can be characterized as an al- 
most rigid shift of oxygen states down in energy by 
about 1 eV for LDA' calculations, while Ni-3d states are 
only slightly modified due to change of Ni-0 hybridiza- 
tion. In other words LDA' calculations lead to the 
change of charge transfer energy A = \sd — e p \ by about 
1 eV. Rather small influence of E^ A on Ni-3d states 
is not surprising, since E^ A for metallic (LDA pro- 
duces metallic state for NiO) electron densities r s =2-6 
are known to be of the order of 1 eV [9] . Further we per- 
form DMFT calculations using LDA and LDA' Hamil- 
tonians, which include all states (without any project- 
ing). DMFT impurity solver used was Hirsh-Fye quan- 
tum Monte-Carlo algorithm [17]. Inverse temperature 
was taken /3 = 5eV _1 (2321 K) and 80 time slices were 
used, with 10 6 Monte Carlo sweeps. The use of very high 
temperature does not lead to any qualitative effects in 
the results, allowing us to avoid unnecessary computa- 
tional efforts. Parameters of Coulomb interaction were 
chosen as typical for NiO [10, 15]: U =8 eV and J=l eV. 
To obtain DMFT(QMC) densities of states at real en- 
ergies, we employed the maximum entropy method [18]. 

In the Fig. 2 we compare the conventional 



LDA+DMFT (upper panel) and consistent 
LDA'+DMFT (lower panel) results for NiO. Dif- 
ferent lines represent partial Ni-3d(^2 3 ) (solid line), 
Ni-3d(e g ) and oxygen 0-2p (dash-dot line) contri- 
butions to density of states. To obtain 0-2p states 
DMFT(QMC) self-energy was analytically continued 
to real frequencies by Pade approximant method. 
For both conventional LDA+DMFT and consistent 
LDA'+DMFT calculations we used H DC of Eq. (5) 
with recalculated on each DMFT iteration step. 
Corresponding values of H DC are 62 eV (n<j=8.7) 
and 58.13 eV (n d =8.2) for conventional LDA+DMFT 
and consistent LDA'+DMFT respectively. The total 
occupancies of Ni-3d states within LDA and LDA' 
calculations were 8.5 and 8.3. 

Within conventional LDA+DMFT we obtain the 
metallic solution, which contradicts experiments. This 
fact can be explained as follows. We already mentioned 
that LDA and LDA' calculations results differ mainly 
by the values of charge transfer energy A = |e d — e p | . In 
fact, we observed [19] that double counting correction 
essentially affects A, or the other way around, the dif- 
ferent values of A require the different values of double 
counting corrections to obtain the same results. In its 
turn, the different values of double counting correction 
can lead either to metallic or insulating solutions for the 
same set of other parameters [10, 19]. 

Once we employ the consistent LDA'+DMFT ap- 
proach, we obtain the charge transfer insulating solution 
for NiO, which agrees well with other LDA+DMFT cal- 
culations for NiO [10, 15] and experiment [20], confirm- 
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Fig. 2. Consistent LDA'+DMFT (lower panel) and 
LDA+DMFT (upper panel) partial densities of states 
for NiO. The Fermi level is at zero energy. 

ing the effectiveness of our approach. Namely, the peak 
at -2 eV which consists almost in equal parts from Ni-3d 
and 0-2p states is nothing else but Zhang-Rice bound 
state (in agreement with Ref. [15]). Lower Hubbard 
band formed mainly from Ni-3d states is located lower in 
energy than Zhang-Rice band. Conducting band is just 
the upper Hubbard band dominated by Ni-3d states. 

As an additional check of consistency of our ap- 
proach we also performed LDA'+DMFT calculations 
for SrVC>3. The results obtained are in good agreement 
with those obtained in Ref. [21], further validating our 
proposed LDA'+DMFT approach as an effective and 
unambiguous method of band structure calculations for 
strongly correlated systems. 

4. CONCLUSION 

In this work we proposed the consistent 
LDA'+DMFT approach, which solves the problem 
of non-uniqueness of the double counting correction. 
By excluding LDA exchange correlation contribution 
for correlated states within the self-consistent LDA 
calculations (e.g. for Ni-3d states) we end up with 
just Hartree like portion of interaction for (strongly) 
correlated states. Then we know exactly, what should 
be subtracted as a double counting correction term, 



while merging LDA' and DMFT. We tested our 
consistent LDA'+DMFT approach, calculating the 
band structure of NiO. We obtained the insulating 
solution without any additional fitting parameters 
and in general agreement with experimental data [20], 
while in other LDA+DMFT works for NiO the double 
counting correction was either treated as an adjustable 
parameter [10], or the special form of double counting 
term was introduced [15] to achieve agreement with 
experiment. 
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n-2-1002) and of the Physics Division of RAS "Strongly 
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